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In  an  interesting  report*  on  rational  approximations  for  the 
l*cx 

function  x  ,  Kopal  suggested  the  following  procedure.  Start  with 
the  approximation 


f«(x)  =  - -  . 

0  1  -  a  +  ax 


Since  the  function  f(x)  =  x  satisfies  a  functional  equation 
f(x)  =  1/  fd/x)  ,  the  function  defined  by  the  equation 


Tq(x)  =  1/$q(1/x) 


is  also  an  approximation  to  x 


1-a 


Now  form  the  new  approximation 


and  repeat  the  entire  process  starting  with  .  It  is  clear  that  we  may 
generate  in  this  manner  a  sequence  of  functions  ,  §2  » • • • 

which  are  rationa?.  As  an  example  of  this  approximation,  we  may  take  x=2 
and  a  =  ^  to  get,  following  Kopal,  the  following  approximations  for  : 


V2) 

?0(2) 

f^(2) 

$3^(2) 

9^(2) 


=  V5  =  1.5335535 
=  5/2  =  1,5000000 
=  17/12  =  1.4666667 
=  24/17  =  1.4117647 
=  577/408  =  1.4142156® 


* 


Zdenek  Kopal,  "The  Approximation  of  Fractional  Powers  by  Rational  Fractions" 
U.  S.  Army  Mathematics  Research  Center,  Report  No.  119*  December  1959. 
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The  formulas  that  one  obtains  by  this  method  are  as  follows  i 


*0^*^  =  1  -  a  +  ox 


♦q(x)  =  o  +  (1  -  a) X 


.  /V  2x  +  o(l  -  a)(l  -  x) 

*1^*^  =  “(i-oy+aox 


¥.  (x)  =  2x(a*x  -  ax) - ^ 

2x  +  a(l  -  a)(l  -  xj 


Kopal  expressed  the  conjecture  that  -•  x'*’  as  n  -•  oo  , 

but  reported  a  proof  only  for  the  case  o  =  .  Our  analysis  below  will 

show  that  this  case  is  the  only  one  in  which  the  conjecture  is  true. 

Kopal  noted  that  when  a  =  the  successive  approximations  to  Jx  ob- 
tained  by  his  method  were  identical  with  those  produced  by  Newton's  method. 

For  the  subsequent  analysis  we  define  the  sequence  of  approximations 
recursively  as  follows : 


«„(x)  = 


1  -  o  +  ox 


¥„(x)  =  1/#  (1/x) 
n  n 


n  =  0,  1,  2,  ... 


♦no./*)  =  (x)  +  ?  (x)]  . 


Fixing  X  and  o  ,  dbfine  X  =  *q(x)  /  t^d/x)  .  We  are  going  to  prove 


that  for  all  n  , 

X  =  ♦  (x)/  §  (1/x)  . 
n  n 

For  n  =  0  this  is  true  by  definition.  If  it  is  true  for  n  ,  then  it 
is  true  for  n  +  1  since 

(x)  i  (x)]  ♦  (x)  +  l/i  (l/x) 

n-t-i  _  n  n  _  n  n 

*  (l/x)  iC#  (l/x)  +  i  (l/x)]  §  (l/x)  +  1/ i  (x) 

n+J.  n  n  n  n 

X  $  (l/x)  +  1/  i  (l/x) 

_  n  n _  _ 

«  (l/x)  +  l/x  *  (l/x) 

n  n 

Now  define  recursively  u^  =  ^q(x)  and  shall 

show  that  u  -  i  (x)  for  all  n  .  When  n  =  0  this  is  true  by  definition 
If  it  is  true  for  n  then  it  is  true  for  n  +  1  since 

n  n 

n 

We  now  see  that  the  sequence  Cu^)  is  the  well-known  Newton-Raphson 

iteration  scheme  for  calculation  of  .  It  converges  if  >  0  and 

\  >  0  From  the  definition  of  u^  ,  we  see  that  these  conditions  are  met 
oc  1  oc 

whenever  x  >  -  and  x  >  - ^  .  Thus  in  these  cases  the  sequence  of 

a  a  —  X 

approximations  #q(x),  ♦^(x),,..  converges  to 
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In  the  case ..  o  =  =  »/x  .  Note  that  any  alteration  in  will 

not  affect  the  proof  of  the  fact  that  “*  •  Only  the  formula 

for  \  will  be  changed. 

Remark  The  functional  equation  mentioned  earlier 

f(x)  =  1/ f(l/x)  , 

which  is  satisfied  by  f(x)  =  x**  ,  has  for  its  general  solution 
f(x)  =  exp  g  (log  x)  where  g  is  any  odd  function  (i.e.,  g(x)  =  -g(-x)) 
This  was  pointed  out  to  me  by  Professor  John  Green,  and  the  proof  is  as 
follows.  Put  X  =  e^  and  h(z)  =  f(e^)  so  that  1  =  f(x)  f(l/x)  = 

=  f(e^)  f(e”*)  =  h(z)  h(-z)  .  Taking  logs,  0  =  logh(z)  +  logh(-z)  , 
and  therefore  logh(z)  is  an  odd  function,  g(z)  .  Thus  f(x)  =  f(e*) 

=  h(z)  =  expg(z)  =  expg(logx)  . 


